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ABSTRACT:  The  attenuation  of  stress  waves  in  samples  of  gabbroic  rock  subjected  to 
shock  loading  up  to  1 1  GPa  is  studied.  We  determine  the  attenuation  coefficients,  ap, 
for  samples  with  different  damage  parameters  under  dynamic  strains  of  2x1  O'7  and  at 
frequencies  around  2  MHz  using  the  ultrasonic  pulse-echo  method.  A  fit  to  the  data 
yields  the  P-wave  spatial  attenuation  coefficient  versus  damage  parameter:  ap  =  40.9  Dp  - 
30.5  Dp2  (db/cm).  Basing  on  O'Connell-Budiansky  theory  a  relation  between  attenuation 
coefficient  and  crack  density  is  derived.  The  predictions  of  ap  and  Q  from  Walsh's 
theory  agrees  well  with  the  experiment  results  for  the  samples  with  different  damage 
deficits  and  the  Q's  we  measure  are  in  the  range  of  10  to  20.  These  very  low  values  give 
rise  to  the  sharp  decrease  in  high  frequency  seismic  energy,  as  the  stress-wave  from  an 
explosion  leaves  the  source  region. 
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OBJECTIVE: 


The  attenuation  of  stress  waves  in  rocks  is  largely  caused  by  internal  friction 
between  crack  surfaces.  When  a  stress- wave  propagates  in  rock  this  induces  differential 
motion  along  the  usually  present  cracks  [Bom,  1941;  Walsh,  1966] .  Generally, 
attenuation  depends  on  stress  amplitude  and  frequency,  as  well  as,  pressure,  temperature 
and  fluid  saturation. 

Experimental  measurements  of  attenuation  of  stress  waves  in  several  rocks  have 
been  extensively  studied  since  the  1940's  in  the  laboratory  using  different  techniques  over 
a  wide-frequency  range  [Bom,  1941;  Toksoz  et  al.,  1979;  Winkler  and  Plona,  1982].  All 
these  studies  concentrated  on  the  relations  between  attenuation  and  parameters  other  than 
crack  density  induced  velocity  deficits.  Because  attenuation  is  directly  related  to  crack 
density,  research  on  influence  of  crack  density  on  attenuation  can  provide  detailed 
information  about  the  mechanism  of  attenuation.  The  present  study  presents  the  first 
experimental  data  of  the  P-wave  attenuation  of  gabbroic  rocks  (San  Marcos,  CA)  which 
are  pre-damaged  by  shock  waves. 

EXPERIMENTAL  TECHNIQUE: 

Sample  Preparation 

The  rock  studied  was  San  Marcos  gabbro  which  has  been  studied  previously 
[Ahrens  and  Rubin,  1993;  Rubin  and  Ahrens,  1991]  The  density  of  San  Marcos  gabbro  is 
2.87  g/cm3,  and  there  is  very  low  initial  crack  density. 

Initially  a  large  gabbro  target  with  dimensions  200x200x150  mm  was  impacted  by 
a  lead  projectile  at  a  velocity  of  1.2  km/s,  the  projectile  had  a  diameter  of  7  mm  and  mass 
of  3  gm.  The  pressure  of  the  shock  wave  in  the  target  was  estimated  to  be  about  1 1  GPa 
at  a  radius  of  0.35  cm  and  about  0.1  GPa  at  radius  of  8  cm  from  the  impact  site  along  the 
center-line  of  the  impact  using  a  power-decay  relation  and  impedance  match  method 
[Ahrens,  1987;  Ahrens  and  Johnson,  1995]. 

The  recovered  target  was  cut  into  1  cm  cubes  with  two  surfaces  being 
perpendicular  to  the  impact  axis.  The  1  cm  cubes  were  polished  to  ±  0.03  mm.  Prior  to 
measurement,  the  samples  were  dried  in  an  oven  under  normal  pressure  at  100°  C  for  24 
hours. 


Damage  Parameter  for  P-wave 

For  the  damage  parameter,  we  follow  the  definition  used  in  Grady  et  al.  [1987] 
and  Ahrens  et  al.  [1995]  which  gives 


\p°  J 


(1) 


where  Dp  is  the  damage  parameter  for  a  P-wave,  Cp  is  the  P-wave  velocity  measured  for 
the  damaged  samples  and  Cjjo  is  the  intrinsic  P-wave  velocity  for  initial  samples. 

Attenuation  Coefficient 

The  ultrasonic  experimental  apparatus  used  in  this  work  is  similar  to  that 
developed  by  Winkler  et  al.  [1982]  for  attenuation  coefficient  measurements  (Fig.  (1)). 
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Fig.  1  Measurement  system 


Suppose  that  the  reflecting  and  transmission  coefficients  of  the  surfaces  between 
the  buffer  and  sample  are  equal  to  that  for  plane- wave  incidence  as  in  Winkler  et  al. 
[1982],  the  changing  amplitude  of  the  wave  is  then  totally  from  intrinsic  attenuation.  Let 
R  be  reflection  coefficient  for  the  interface  between  the  coupling  buffer  and  sample  and 
L/2  be  sample  thickness,  A(f)  and  B(f)  be  the  frequency-dependent  amplitudes  of  the 
pulses  reflected  from  the  surface  A  and  B  (in  Fig.  1)  of  the  sample,  respectively.  The 
attenuation  coefficient  obtained  from  the  two  reflected  stress  waves  is  expressed  as 
[Winkler  and  Plona,  1982] 


a(f)  = 


8.686, 

- In 

L 


A(f) 

B(f) 


(2) 


where  the  unit  of  a(f)  is  db/cm  when  the  unit  of  L  is  cm  and  R  is  expressed  as 


R  =  CpP  CpcPc 

CpP  +  CpcPc 


(3) 


where  Cp  and  p  are  the  P-wave  velocity  and  the  density  of  the  gabbro  sample  measured, 
respectively.  Cpc  and  pc  are  the  P-wave  velocity  and  density  of  the  coupling  buffer 
(Lucite).  From  ultrasonic  measurements,  Cpc  is  2.68  km/s  and  pc  is  1.19  g/cm3. 

A  piezoelectric  transducer  (Panametrics,  Model  102)  is  used  as  the  pulse 
generator  and  the  receiver.  The  transducer's  driver  is  a  Panametrics  505UA 
pulser/receiver.  The  recorded  signals  are  used  to  calculate  the  magnitude  of  each 
frequency  component  of  each  pulse  using  a  Fast  Fourier  Transformation  (FFT). 

PRELIMINARY  RESEARCH  RESULTS  AND  ANALYSTS: 

A  typical  signal  recorded  for  the  attenuation  measurement  is  shown  in  Fig.  2,  in 
which  the  first  signal  is  reflected  from  the  surface  A  of  the  sample  and  the  second  from 
the  surface  B,  and  Ti  and  T2  are  the  data  length  used  in  FFT.  We  also  measured  the 
strain  in  the  rock  produced  by  the  ultrasonic  transducer.  The  strain  level  is  10'7  and, 
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the  present  data  are  taken  in  the  ultrasonic  region,  the  results  are  also  applicable  to 
seismic  frequencies. 
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Fig.  2  Typical  ultrasonic  record 

A  typical  result  of  spectral  analysis  is  shown  in  Fig.  3,  it  is  found  that  most  of 
energy  of  the  stress  wave  generated  by  the  transducer  is  concentrated  in  the  frequency 
range  between  1.5  and  3.5  MHz.  Using  Eq.  (2),  the  attenuation  coefficients  for  the 
samples  with  differential  damage  parameters  have  been  evaluated  as  shown  in  Figs.  (4) 
(5)  and  (6).  From  Fig.  (4),  we  can  see  that  the  attenuation  coefficients  for  samples  with 
different  damage  parameters  increase  approximately  linearly  with  frequency.  From  Fig. 
(5),  the  relation  between  attenuation  coefficient  and  damage  parameters  for  2  MHz  are 
fitted  as 

otp(db/cm)  =  40.9Dp~30.5Dp2  (4) 
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Fig.  3  Typical  spectral  amplitude  of  signals 
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Fig.  4  Experiment  results  of  attenuation  coefficient 
with  frequency  for  samples  with  different 
damage  parameters 
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Fig.  5  Relation  between  damage  deficits  and  Fig.  6  Quality  factor  versus  damage  deficits 
attenuation  coefficient 


O'Connell  and  Budiansky  [1974]  established  the  relation  between  velocity  and 
crack  density  follows  as 


lCroJ 


(l-vXl-HVo)f  16(1 -vV 
(l  +  v)(l-v0)[  9(1  - 2v)  J 


(5) 

(6) 


where  v  and  Vo  are  effective  and  intrinsic  Poisson’s  ratio,  respectively.  £  is  crack  density 
which  is  defined  as 

£=  N  <  a3  >  (7) 

where  N  is  the  number  of  cracks  per  unit  volume  and  a  is  the  half-length  of  cracks. 

From  the  definition  of  damage  parameter,  Eq.  (1),  and  the  relation  between 
velocity  and  crack  density,  Eqs.  (5),  (6)  and  (7),  the  relation  between  Dp  and  e  is 
approximately 


Dp  =  2.4e- 1.2e2  (8) 

From  Eqs.  (4)  and  (8),  an  approximate  expression  relating  exp  to  £  is 

ap  =  96.1e(1-2.2e  +  1.69e2-0.4e3)  (9) 

Figure  7  presents  experimental  results  of  attenuation  versus  crack  density  and  the 
fitted  curve.  It  can  seem  that  increasing  rate  of  attenuation  coefficient  versus  crack 
density  decreases  with  crack  density.  This  can  be  explained  using  the  relations  among 
attenuation  coefficient,  crack  density  and  half-length.  From  the  definitions  of  the 
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attenuation  coefficient  and  crack  density,  the  attenuation  coefficient  is  related  to  the  area 
of  crack  surface,  this  means  ap  «  a2,  therefore,  from  the  definition  of  crack  density,  we 
have 


(10) 


where  b  is  a  constant.  Therefore,  the  ratio,  — ,  must  decrease  with  increasing  of  crack 

e 

density.  This  is  shown  by  the  experimental  results  of  Fig.  7.  Moreover,  in  Fig.  6  we 
suppose  b  is  a  constant,  the  average  half-length  of  cracks  can  be  estimated  as  shown  in 
Fig.  8. 


e 

Fig.  7  Relation  between  attenuation 
coefficient  and  crack  density 


e 


Fig.  8  Estimated  average  half-length  of 

cracks  versus  crack  density  parameter 


Basing  on  the  concept  that  frictional  dissipation  at  crack  surfaces  in-contact  slide 
relative  to  one  another  during  the  passing  of  stress  wave,  Walsh  [1966]  developed  an 
expression  for  the  attenuation  coefficient  that  can  approximately  be  expressed  as 


a  =Hfe 


K(l-v)(l-2v) 

K0(l-2v)(l-2v„) 


(11) 


where  H  is  considered  to  be  a  constant  and  f  is  frequency.  Only  one  parameter,  H,  needs 
to  be  determined.  If  let  e  =  0.1,  from  Eqs.  (6),  (8),  the  relation  between  K  and  Ko 
[O'Connell  and  Budiansky,  1974]  and  experimental  results  of  attenuation  coefficient.  H  is 
evaluated  to  be  60.7db/cm  for  f  =  2  MHz. 

The  calculated  results  from  Eq.  (1 1)  is  shown  in  Fig.  7.  It  is  obvious  that  the  form 
of  Walsh's  theory  is  in  good  agreement  with  the  experimental  results  and  the  relation 
between  ap  and  e  is  not  linear.  This  result  demonstrates  that  Walsh's  theory  can  be  used 
to  evaluate  attenuation  coefficient  in  rocks  with  high  crack  densities  such  as  occur  in  the 
very  close  vicinity  of  contained  explosions. 
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CONCLUSIONS  AND  RECOMMENDATIONS: 


Some  20  samples  of  San  Marcos  gabbro  cut  from  a  shock  loaded  target  are  used 
to  measure  attenuation  coefficients  for  different  damage  levels  experimentally  using  the 
ultrasonic  method.  We  used  O'Connell  and  Budiansky's  theory  for  relating  crack  density 
to  elastic  constant  with  various  damage  parameters.  Walsh’s  theory  [1966]  was  used  to 
predict  the  attenuation  coefficients  of  the  samples  with  different  crack  density.  Our  main 
conclusions  are: 

(1)  The  P-wave  attenuation  coefficient  and  the  damage  parameters  are  given  by  Eq.  (4). 

(2)  Based  on  the  O'Connell  and  Budiansky's  theory,  die  relation  between  the  crack 
density  and  attenuation  coefficient  can  be  expressed  approximately  as  Eq.  (9). 

(3)  Basing  on  Walsh’s  theory,  we  evaluated  the  attenuation  coefficient.  The  predictions 
of  attenuation  coefficient  fit  experimental  results  as  shown  in  Fig.  (7).  Moreover,  the  Q’s 
of  damaged  rocks  are  low,  10-20.  This  gives  rise  to  the  very  strong  attenuation  of  the 
high  frequency  portions  of  the  seismic  signals  from  contained  explosions. 
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